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(Dated: September 6, 2008) 
Abstract 

Using samples of (5.93 ± 0.10) x 10^ T(3S') decays and (9.11 ± 0.14) x 10^ T(2S') de- 
cays collected with the CLEO detector, we report improved measurements of the branching 
fractions for the following five transitions: ^(T(35)^T(lS')7r+7r-) = (4.46 ± 0.01 ± 0.13)%, 
H(T(2S)^T(lS)7r+7r-) = (18.02 ± 0.02 ± 0.61)%, B{T{3S)^T{lS)'!T^7r°) = (2.24 ± 0.09 ± 0.11)%, 
H(T(2S)^T(lS)7rV) = (8.43±0.16±0.42)% and B(T(35)^T(2S)7rV) = (1.82±0.09±0.12)%. 
In each case the first uncertainty reported is statistical, while the second is systematic. 
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Hadronic transitions among heavy quarkonium states provide an excellent testing ground 
for non-perturbative Quantum Chromodynamics (QCD) They are generally understood 
to proceed by the emission and hadronization of low momentum gluons [2,], and their inves- 
tigation is one of few possible laboratories for the study of the low-g^ hadronization process. 
The study of such transitions in the bottomonium (bb) system is particularly advantageous 
because of the non-relativistic nature of the system and the richness of the spectrum of 
states below open-bottom threshold. (See Figure [H) 

For the first 22 years after the observation of hadronic transitions among bottomo- 
nium states by LENA [3] and CUSB [4], only the six tctt transitions among the vector 
T{nS) bottomonia were known. CLEO has recently observed three other examples of 
hadronic transitions in bottomonium: X6i,2(2P)— >ct;T(lS') x;,(2P)^7r7rXft(lP) \^ and 
T {2S)—>-r]T {IS) Very recently, the BaBar Collaboration has reported new measure- 
ments of several hadronic transitions in the bottomonium system using bottomonium states 
produced in ISR while running at the T{4S) resonance [8|. 

In this Article we report improved measurements of the branching fractions for nn transi- 
tions among the vector states of the bottomonium system. Dipion transitions from T(35') to 
the lower vector states (T{2S), T^ISY) and from T(2S') to T(IS') have been of interest ever 
since their first observation in 1982 [9|. There has recently been a resurgence of interest in 
dipion transitions following the observation by Belle [ll[ and BaBar [l2| of dipion transitions 
from the bottomonium resonances T{4S) and T{5S) to lower bottomonium states, of the 
new state Y(4260) ll3fL_to J/ip, and also the observation by BES and CLEO of similar tran- 
sitions of V^(3770) [ij, ll5| . Additional motivation to update measurements of the branching 
fractions for dipion transitions among bottomonium states below open-bottom threshold is 
presented by the prospects of using T(35'), T(2S')— s>7r7rT(lS') as a clean source of tagged 
T(IS') to study exclusive T(15') decays, including searches for invisible decay modes [l6j. 

A well-known feature of the tttc transitions in bottomonium is that the invariant mass of 
the dipion system in T(35')^7r7rT(15') differs greatly from that produced in other known 
dipion transitions in bottomonium and in charmonium P]. Theoretical interest in these 
invariant mass distributions has been substantial, and several attempts to describe them 



have been made since the first observation of dipion transitions [10|. Analysis of the dipion 
invariant mass shapes in transitions from bottomonium states above open-bottom threshold 
show similar interesting features 11, T^]. A detailed analysis of the dipion invariant mass 



shapes including the extraction of the matrix elements for the transitions considered in this 
Article was performed using the most recent CLEO data and appears in Ref. [13]. The 
results of that matrix element determination are used in the present analysis to properly 
determine the detection efficiency. 

For this work, data samples were collected using the CLEO III f[§\ detector at the Cornell 
Electron Storage Ring. These samples included (5.93 ± 0.10) x 10^ T(3S') collected on the 
T(3^) resonance, at = 10.355 GeV and (9.11 ± 0.14) x 10^ T(2^) decays collected on 
the T(25') resonance, at ^/s = 10.025 GeV. Charged particle tracking is done by a 47- 
layer drift chamber and a four-layer silicon tracker immersed in a 1.5 T solenoidal magnetic 
field. Photons are detected using an electromagnetic calorimeter consisting of 7784 CsI(Tl) 
crystals in a projective barrel geometry. 

In this analysis, we study the transitions both inclusively (in which case we detect only 
the pair of charged pions) and exclusively (in which case we detect, in addition to the 
charged or neutral pair of pions, the decay of the daughter T{nS) state to either f^~^fj,~ or 
e~^e~). In each case, the primary quantity used to identify our observation of the dipion 
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FIG. 1: The spectrum of bottomonium states below BB threshold for different J^*^ combinations. 
Established bottomonium states are indicated by the solid horizontal lines, while those indicated 
by dashed horizontal lines have never been observed. Arrows connecting various bottomonium 
states represent hadronic transitions that have been observed. 



transitions of interest is the mass recoiling against the dipion system. This may be most 

)2, where Erm is the 



simply defined in terms of the formula Mrecoii 



energy in the center-of-mass system, ^^tttt and Ptttt are the energy and three-momentum of 
the dipion system, respectively. For a dipion system produced in the transition from T{nS) 
to T(IS') (T(25')), Mrecoii will be equal, within detector resolution, to the T(IS') (T(2S')) 
mass, 9.460 (10.023) GeV/c^ [19]. Randomly selected pion pairs from hadronic events do 
not peak at all, but form a smooth combinatoric background. 

In the inclusive analysis, we select two charged tracks that originate within 5 cm in the 
beam direction and 5 mm in the transverse direction from the center of the interaction region. 
The vertex requirements drastically reduce the likelihood that the vr candidate tracks were 
produced in interactions between the e"*" or e~ beams and the beampipe, residual gas, or 
other material. Monte Carlo studies show that low momentum kaon or muon pairs arising 
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from generic T{nS) or QED processes do not produce peaking backgrounds in the recoil mass 
spectrum. Therefore, we did not require positive identification of these tracks as pions. The 
mass recoihng against the two tracks, assumed to be pions, was required to be greater than 
9.0 GeV/c^ 

In the exclusive analysis we require events to have two high momentum tracks with an 
invariant mass of 9.2 - 9.7 GeV/c^, consistent with M(T(15)), or greater than 9.9 GeV/c^, 
consistent with M{T{2S)). These tracks must have an angle with respect to the beam 
direction, 9i, which satisfies | cos^^j < 0.82, a region in which the acceptance is relatively 
uniform. We apply no further track quality criteria for the lepton candidates, and we do not 
attempt to distinguish to which dilepton final state (electron or muon) the T(15') candidate 
has decayed. Since in the exclusive analysis we reconstruct the full event, we require the sum 
of energies of all final state particles to be greater than ^/s — 200 MeV(— 300 MeV)) for the 
T(35') (T(25')) analyses. The more stringent requirement on the T(35') energy conservation 
was necessary to remove possible contamination in the signal sample due to cascades through 
the T{2S). 

In the charged exclusive case, we require that in addition to the dilepton candidate, 
events have a pair of low momentum tracks that satisfy the same requirements as the pion 
candidates in the inclusive analysis. The dilepton invariant mass requirement alone provides 
a nearly background- free sample, and imposition of additional criteria in order to identify the 
tracks as leptons only leads to larger systematic uncertainties and reduced signal efficiency 
without much improvement in signal quality. An additional requirement is imposed to 
remove radiative Bhabha events (e+e"— >e+e~7) in which the 7 converts in the inner material 
of the detector or the beampipe, producing an e'^e~ pair that can fake the transition vr+vr". 
In such events a small angle between momentum of the conversion pair and one of the two 
high-momentum leptons is favored. Essentially all of this background is removed by the 
requirement that this angle be greater than 0.15 radians. Finally, in addition to the four 
tracks (two candidates, two candidates) we allow events to have one additional track, 
which prevents the loss of otherwise good events due to failures in pattern recognition or to 
other spurious track candidates. Monte Carlo (MC) studies show that this allowance does 
not contribute any peaking background in the region of interest. 

In the neutral dipion analysis, we require that in addition to the dilepton candidate, events 
contain four or five showers in the calorimeter. Each of these showers must have an energy of 
at least 50 MeV and have angles relative to the beam axis such that | cos(^^^)| < 0.804, where 
7 reconstruction is best. None of these showers may be matched to either charged track 
in the event. Showers satisfying these criteria are then paired to produce 7r° candidates. 
Pairs with invariant masses within 50 MeV/c^ of the nominal 7r° mass must further satisfy 
the requirement that the mass pull - the normalized deviation from the nominal mass, 
(M^^ - M^o)/a^^ - be in the range (-4.0,2.5) for T(n5)^T(15)7r°7r° and (-5.0,3.0) for 
T(35')— >7r°7r°T(25'). The pair-mass resolution, cr^^, which is calculated event by event 
based on the calorimeter energy deposits, is typically between 5 and 7 MeV/c^. The 
candidates satisfying this condition are then subject to a mass-constrained fit in order to 
improve the 7r°7r° recoil mass resolution. 

We have also allowed one additional spurious shower to be present in the calorimeter, in 
order that we do not needlessly remove events due to spurious signals in the calorimeter. 
Because we allow up to five showers in the calorimeter it is possible for an event to have 
more than two combinations of showers that satisfy the 7r° candidate requirements. In such 
cases, the combination of 7r° pairings having the smallest sum of squared mass pulls. 
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In order to evaluate detector acceptances and efficiencies and to study backgrounds to 
the signal processes, MC samples were generated for several different event types. Generic 
T(3S') and T{2S) decays and continuum processes (such as e"'"e~— >-rr) at center of mass 
energies equal to the masses of the two states were simulated using the routine QQ MC 
samples were also generated for the signal dipion transitions and for individual background 
channels using EvtGen 21|. Each sample was then passed through a GEANT-based 13 



detector simulation. Generic and continuum MC samples contained approximately one and 
five times the actual integrated luminosity taken at each of the resonances, respectively. 

For the study of acceptance and efficiency, separate signal MC samples for inclusive and 
exclusive analyses were created. In order to take advantage of the matrix element analysis 
previously performed by CLEG [13], the signal MC samples were generated according to 
phase space and then weighted according to the square of the matrix elements. In the 
exclusive analyses, a further weighting factor of 1 + cos^ (where 6^ is the lepton angle 
relative to the beamline in the rest frame of the daughter T) was applied. This assumes 
a negligible D-wave component in the vrvr transition. For the inclusive analyses, 200,000 
events were generated for each of T{3S) and T(25') decaying by 7r+7r~T(15'), where T(IS') 
was decayed generically. For the exclusive analyses, 500,000 events were generated for each 
of the five transitions, with the daughter T(IS') or T{2S) decaying equally to e~^e~ and 

In all analyses, the distribution of the mass recoiling against the vrvr pair for accepted 
data and MC events is used to evaluate signal yield and efficiency, respectively. Recoil mass 
distributions for all the transitions observed in our data are shown in Figs. [2] and [31 

For the inclusive analysis, there is a large combinatoric background because the analysis 
involves only combining pairs of charged pions, which are prolifically produced in T decays. 
The background due to these is smooth and has been fitted to a third-order polynomial, 
and then subtracted in order to evaluate the yield for the signal process of interest. A study 
of continuum data taken below the T{3S) and T(2S') resonances, and both continuum and 
generic T{3S) and T(2S') MC simulations (which include decays of daughter T(IS') states) 
reveal no peaking background in tt+tt" recoil mass. 

For the exclusive charged analysis, background events arise either from non-vrvr hadronic 
transitions to the daughter T state in which the hadrons involved in the transition fake the 
signal due to poor reconstruction or noise in the detectors, or from udsc quark pair pro- 
duction from the continuum. Backgrounds arising from these process should be negligible, 
because of the lack of significant non-vrvr hadronic transitions in the bottomonium system, 
and because of the required detection of a high mass dilepton. Continuum data taken be- 
low T{3S) and T(2S') resonances and both continuum and generic MC simulations (with 
simulated vr+vr^ transitions removed) have been analyzed, and this expectation is con&med. 

For the neutral analysis, there are also a number of possible background processes 
to consider. One such process is the decay chain T(3S')— >-7Xf)(2P), X6(2-P)— ^7^(25'), 
T(2S)^7Xfe(lP), Xb{lP)^lT{lS), T(l^)^£+r . The resulting 47+£+r final state has an 
overall branching fraction of approximately 0.2%, while the branching fractions of the signal 
processes are of order several percent [l9| . Such a process is furthermore unlikely to produce 
a fake signal because the four photons produced in the cascade would need to combine to 
produce two good 7r° candidates. Similarly, the two-photon cascades T(35')^7X;)(2P, IP), 
followed by XbC^P, lP)^-fT{lS) and T{1 S)^i+i~ could potentially fake the vr^vr" signal 
with the addition of two spurious clusters in the calorimeter. To evaluate these small but 
possibly pernicious backgrounds, special MC samples, each of 50,000 events, were produced 
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FIG. 2: vr^vr" recoil mass distributions for the four charged dipion transition analyses, (a) 
T{3S)^TT+TT-T{IS) inclusive, (b) T(2S)^7r+7r-T(15) inclusive, (c) T(35)^7r+7r-T(15) exclu- 
sive and (d) T(25')^7r+7r~T(15) exclusive. The data are represented by symbols with uncertainty, 
while the histograms overlaid in each inclusive plot represents background-subtracted data. Note 
the logarithmic scale for the exclusive analyses. 



for both the four and two photon cascade processes described above. Analysis of these 
samples show that such processes contribute negligible background. Finally, transitions in- 
volving 7] could as well, but because of the tiny branching fraction for t] transitions [2], they 
have been assumed to produce no background. 

To determine the branching fractions for the charged exclusive analyses, we used a cut- 
and-count method in which the background underneath the peak in the recoil mass distribu- 
tion was subtracted by means of 20 MeV wide sidebands separated by 15 MeV from the 40 
MeV wide signal region. We analyzed the inclusive recoil mass distributions similarly, first 
subtracting the smooth combinatoric background followed by cut-and-count using the same 
signal and sidebands as in the exclusive analysis. This method minimizes the systematic un- 
certainties associated with the MC signal shape, and corrects for any inaccuracies resulting 
from the use of the background function in the vicinity of the signal peak. Finally, from the 
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FIG. 3: vr'^vr" recoil mass distributions in the neutral dipion transition analyses, (a) 
T{3S)^tt\°T{1S), (b) T{3S)^tt\^T{2S) and (c) T{2S)^7:^7r^T{lS). The data are represented 
by the symbols with errors, while the histograms represent the result of the fit using the MC shape 
and a linear background function. Note the logarithmic scale. 



resulting data recoil mass histograms, we obtain the data yield, and from the signal MC, we 
obtain the efficiency by dividing the MC yield by the amount of MC generated. 

For the neutral analyses we instead obtained the efficiency-corrected data yield directly 
by fitting the data recoil mass histogram using the recoil mass histogram shape found from 
neutral signal MC samples, allowing for a first degree polynomial background contribution. 

From the data yields and MC efficiencies, we then calculate the branching fraction, using 



B{T{nS)^T{mS)7rn) 



in the case of the exclusive analyses and 



A^cxc/« 



[N{T{nS))2B{T{mS)^fi+l2-)] 



(1) 



B{T{nS)-^T{lS)7r+7T- 



^ ' mc/ 



(2) 



in the case of inclusive studies. 

The branching fractions obtained are summarized in Tables HI and HTl In order to evaluate 
the exclusively-measured branching fraction, we have assumed lepton universality. We have 
therefore used twice the PDG average value of i3(T(lS')— = 2.48±0.05% to normalize 
our results for the transitions terminatijig in T(IS'), and twice the recent CLEO measurement 
of i3(T(2^)^/i+/i-) = 2.03 ± 0.08% i23j for the one terminating in T(25). 

Systematic error contributions are summarized in Table IIIII For event reconstruction in 
the exclusive analyses, a systematic uncertainty of 1.2% per pair of charged pions and 3.2% 
per pair of neutral pions was assessed. These systematic uncertainties were evaluated by 
comparing the ratio of event yield in the standard exclusive analysis to the yield obtained 
using an analysis that depends on the reconstruction of only one vr^ or 7r°. From this ratio 
a per-vr^ or 7r° uncertainty was obtained, and doubled to give the relative uncertainty for 
finding the pair. For the exclusive analyses, a systematic uncertainty of 1.0% per lepton pair 
was similarly obtained. For the inclusive analyses, based on tracking studies in a variety 



8 



TABLE I: Results of the branching fraction measurements for charged dipion transitions 
T(3S')^7r^7r^T(lS') and T(2S')^7r+7r^T(lS'). The first uncertainty in the branching fraction 
is the statistical uncertainty, while the second is systematic. The averages listed are weighted aver- 
ages of the inclusive and exclusive measurements which take into account correlation of systematic 
uncertainties between them. 



Analysis 


Data Yield 


Efficiency (%) 


B (%) 


3S Excl. 


5215 ± 72 


39.7 ±0.1 


4.47 ±0.06 ±0.18 


3S Inch 


184760 ± 430 


69.9 ±0.2 


4.46 ±0.01 ±0.14 


Average 






4.46 ±0.01 ±0.13 


2S Excl. 


26417 ± 163 


32.0 ±0.1 


18.26 ±0.11 ±0.81 


2S Inch 


824418 ± 908 


50.3 ±0.1 


17.99 ±0.02 ±0.59 



Average 18.02 ± 0.02 ± 0.61 



TABLE II: Results of measurements of the branching fractions for neutral dipion transitions. The 
middle column lists the efficiency-corrected data yield, which is obtained as described in the text. 
The statistical uncertainty presented accounts for both the data and finite MC statistics, and the 
second uncertainty reflects the remaining systematic contributions. 



Analysis 


Efficiency-corrected Yield 


Bi%) 




6584 ± 274 


2.24 ±0.09 ±0.11 


3S^257r07rO 


4391 ± 207 


1.82 ±0.09 ±0.12 




38069 ± 727 


8.43 ±0.16 ±0.42 



of neutral and charged multiplicity environments, we conservatively assign a systematic 
uncertainty of 2.4% per pair of charged pions. 

For the T(35') analyses, a common relative uncertainty of 1.7% due to the uncertainty in 
the number of T(35') produced; for T{2S), the corresponding uncertainty was 1.5%. This 
class of systematic uncertainty derives primarily from the uncertainty in the knowledge of 
the integrated luminosity accumulated at each of the resonances. 

Three sources of systematic uncertainty produce relatively large contributions in some 
cases. Uncertainties due to modelling of the dipion dynamics were studied by varying the 



MC weighting according to the uncertainties of the matrix elements reported in Ref. [17 
and by studying the resulting reproduction of the dipion invariant mass. In the case of the 
T(35')— s>T(2S')7r°7r° and both analyses of T(2S')^T(lS')7r+7r~ estimated systematic contri- 
butions of 1.4 — 2.3% were obtained. For the other transitions, the systematic uncertainty 
due to modelling was much smaller. 

The exclusive samples were all divided into /i"*"/!" and e~^e~ subsamples for the purpose 
of studying the difference between reconstruction of these two leptonic channels. From the 
difference in branching fractions obtained from the two subsamples, a relative systematic 
uncertainty of 2.5% due to lepton type was obtained. The exclusive analyses carry an 
additional systematic uncertainty from the branching fractions for the decays of T(IS') and 
T{2S) to dileptons. 
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TABLE III: Summary of relative systematic uncertainties on the measurement of the branching 
fractions, expressed in percent. Systematic uncertainties for the inclusive and exclusive analyses 
of charged dipion transitions are separately hsted in the table. 







T(35)- 








T{2S)^ 




Contribution 


T(15)7r+^- 


T(15)7r+7r- T{1S)tt^tt^ 


T(25)7r07r0 


T(lS)7r+7r- 


T(15)^+^~ 


T(15)7r0^o 




Excl. 


Incl. 






Excl. 


Incl. 




7r±/7rO 


1.2 


2.4 


3.2 


3.2 


1.2 


2.4 


3.2 


£ Tracks 


1.0 


N/A 


1.0 


1.0 


1.0 


N/A 


1.0 


Luminosity 


1.7 


1.7 


1.7 


1.7 


1.5 


1.5 


1.5 


£ Type 


2.5 


N/A 


2.5 


2.5 


2.5 


N/A 


2.5 


MC Modelhng 


0.2 


0.4 


0.5 


2.2 


2.3 


1.4 


0.2 


U BR 


2.0 


N/A 


2.0 


4.2 


2.0 


N/A 


2.0 


Other Sources 


0.35 


0.8 


1.0 


1.0 


0.1 


0.8 


1.0 


Total 


4.0 


3.1 


5.1 


6.6 


4.5 


3.3 


5.0 



Uncertainties due to the choice of analysis requirements, MC statistics, side band range 
choices, etc., were all much smaller in each case compared to the other systematic uncer- 
tainties. The overall relative systematic uncertainty is obtained by adding all contributions 
in quadrature. The complete array of systematic uncertainties for all seven analyses appears 
in Table [TTll 

We have measured the charged dipion branching fractions both inclusively and exclusively. 
Our final result for the measurement of each of these branching fractions is a weighted average 
of the two independent results, which we have caluclated using a toy Monte Carlo method [2J] 
in which we have properly accounted for the correlation between the various contributions 
to the systematic uncertainties that are applied to the two analyses. That is, the luminosity 
uncertainties on the exclusive and inclusive results are fully correlated, while the statistical 
uncertainties are uncorrelated, as are any uncertainties unique to either the exclusive or the 
inclusive analysis. These average values are i3(T(3S')— i>T(lS')7r+7r~) = (4.46 ±0.01 ±0.13)% 
and S(T(2^)^T(15)7r+7r-) = (18.02 ± 0.02 ± 0.61)%. 

It is interesting to compare the branching fractions for the vr^vr^ transitions to those for 
the corresponding vr+vr" branching fractions that we have measured. Isospin conservation 
requires that the square of the matrix elements for the 7r°7r° transitions be half that of 
the TT+vr^ transitions. Phase space for the two types of transitions also differs slightly and 
modifies this expectation, such that the expected ratio i3(7r°7r°T(15'))/i3(7r+7r^T(15')) is 
0.53 for the transitions from T(2S') and 0.51 for transitions from T(3S'). Combining our 
neutral and charged results, and taking into proper account correlations and cancellations 
among individual systematic errors, we obtain ratios of 0.462 ± 0.037 for T(2S') transitions 
and 0.501 ± 0.043 for the transitions from T(3S'). 

In summary, we have reported improved measurements of five of the six dipion transitions 
among the lower-lying bottomonium vector states T(3S'), T(2S') and T(15'). Each of the 
measurements is more precise than those made by any previous experiment, and also more 
precise than the current PDG world average [l9| . 
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